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We analyze the end-point region of the photon spectrum in semi-inclusive radiative decays of 
very heavy quarkonium (ma 2 3> Aqcd)- We discuss the interplay of the scales arising in the 
Soft- Collinear Effective Theory, m, m(l — z) 1 ^ 2 and m(l — z) for z close to 1, with the scales of 
heavy quarkonium systems in the weak coupling regime, m, ma B and ma 2 . For 1 — z ~ a 2 only 
collinear and (ultra)soft modes are seen to be relevant, but the recently discovered soft-collinear 
modes show up for 1 — z <C a 2 . The S- and P-wave octet shape functions are calculated. When 
they are included in the analysis of the photon spectrum of the T(IS') system, the agreement with 
data in the end-point region becomes excellent. The NRQCD matrix elements (l 3 5i|08( 1 <S'o)|l 3 £'i) 
^2 . and (l 3 S , i|0 8 ( 3 Pj)|l 3 5'i) are also obtained. 

o . 

t PACS numbers: 13.20.Gd, 12.38.Cy, 12.39.St 

I—* 

'. I. INTRODUCTION 

Effective field theories (EFT) have proved extremely useful in the field of strong interactions. Applications to high 
energy processes in QCD involving very energetic partons 0, however, have been elusive until recently. Important 
C<") ■ features of a suitable EFT for such processes were outlined in 0, which led to the development of the so called 
' Soft-Collinear Effective Theory (SCET) 0[IH3 (see for a pedagogical introduction). 

SCET has generated high expectations. Indeed factorization proofs appear to greatly simplify and power corrections 
seem to become under control. In addition, a large number of potential applications is envisaged J). Among those, 
exclusive and semi-inclusive i?-decays have deserved special attention due to the necessity to have a good control on 
the hadronic effects in order to extract the CKM matrix elements from the abundant -B-factory data. 

SCET was originally formulated in terms of soft, collinear and usoft modes. Later, it was realized that two possible 
scalings for collinear modes were relevant and the terminology SCETj and SCETn was introduced. Recently a new 
mode, called soft-collinear, has been claimed to be necessary Q (see Q for the latest discussions). It is often assumed 
<D , that some of the momentum components of these modes have typical sizes ~ Aqcd or even smaller. 

One of the difficulties that one faces in i?-physics is that the bound state dynamics of the initial JS-meson is 
^ , dominated by the scale Aqcd, and hence a weak coupling analysis is not reliable. Therefore, the interplay of the 
initial bound state dynamics with final state modes of momentum components of the order Aqcd (or smaller) is 
difficult to figure out. We advocate here that a very heavy quarkonium in the initial state may provide an excellent 
theoretical tool to shed light on this issue, since the bound state dynamics occurs at weak coupling and it is amenable 
of a detailed analysis. We shall illustrate this point by analyzing the end-point region of the photon spectrum in 
inclusive decays of very heavy quarkonium. 

Semi-inclusive radiative decays for the T(IS') have already been discussed in the framework of SCET |sL ITot fTT| . 
SCET has been used to put forward factorization formulas and to resum Sudakov logarithms. An improved description 
of data ^3 with respect to earlier approaches 01 has been achieved. However, the bound state dynamics, which is 
relevant for the evaluation of the octet shape functions, has not been studied in detail, but rather modeled by analogy 
with B-meson systems 14], which is a doubtful approximation. We shall calculate here the octet shape functions 
under the assumption that the bottom quark is sufficiently heavy as to consider T(15') a Coulombic state. This 
assumption appears to be self-consistent in the calculation of the spectrum [151 Il6l Il7l Il8| , and decay and production 
currents • We observe that the factorization scale dependence of the shape functions is sensible to the bound state 
dynamics and discuss its cancellation. 

We distribute the paper as follows. In section [H] we calculate the photon spectrum at the end-point region in 
the weak coupling regime. We do so by matching first QCD to Non-Relativistic QCD (NRQCD) |2(j +SCETi, next 
NRQCD+SCETi to Potential NRQCD (pNRQCD) H3+SCET n , and finally carrying out the calculations in the later 
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EFT. We confirm the factorization formulas [ToL Hl| and obtain the octet shape functions. In section HTT1 we apply 
our results to the T(15) system and obtain a very good description of the experimental data ^2] m the end-point 
region. In section ITvl we discuss the interplay of the several scales in the problem, in particular the emergence of a 
soft-collinear mode. Section is devoted to the conclusions. In the Appendix we present results for NRQCD octet 
matrix elements of the l 3 5i state, which follow from those of section ITT1 

II. THE END-POINT REGION OF THE PHOTON SPECTRUM 



We start from the formula given in 




^ = z^LlmT(z) T(z) = -if fxe-**-* (n 2S+1 L,j\T{ J^x) J v (0)}\n 2S+1 Lj) rff (1) 

where J^x) is the electromagnetic current for heavy quarks in QCD and we have used spectroscopic notation for 
the heavy quarkonium states. The formula above holds for states fulfilling relativistic normalization. In the case 
that non-relativistic normalization is used, as we shall do below, the rhs of either the first or second formulas in (1) 
must be multiplied by 2M, M being the mass of the heavy quarkonium state. At the end-point region the photon 
momentum (in light cone coordinates) in the rest frame of the heavy quarkonium is q = (<7+, g_, q±) = (zM/2, 0,0) 
with z ~ 1 (My/1 — z <C M) . This together with the fact that the heavy quarkonium is a non-relativistic system fixes 
the relevant kinematic situation. It is precisely in this situation when the standard NRQCD factorization (operator 
product expansion) breaks down p2^| . The quark (antiquark) momentum in the QQ rest frame can be written as 
V = (Po 5 p)- Po = 77i + Zo, p = 1; Zo, 1 <C to, to being the mass of the heavy quark (M ~ 2m). Momentum conservation 
implies that if a few gluons are produced in the short distance annihilation process at least one of them has momentum 
r = (r + , r_ ,r±), r_ ~ M/2 ; r + , r± <C M, which we will call collinear. At short distances, the emission of hard gluons 
is penalized by a s (m) and the emission of softer ones by powers of soft scale over M. Hence, the leading contribution 
at short distances consists of the emission of a single collinear gluon. This implies that the QQ pair must be in a 
color octet configuration, which means that the full process will have an extra long distance suppression related to 
the emission of (ultra)soft gluons. The next-to-leading contribution at short distances already allows for a singlet QQ 
configuration. Hence, the relative weight of color-singlet and color-octet configurations depends not only on z but also 
on the bound state dynamics, and it is difficult to establish a priori. In order to do sc^ it is advisable to implement 
the constraints above by introducing suitable EFTs. In a first stage we need NRQCD |2(j, which factors out the scale 
to in the QQ system, supplemented with collinear gluons, namely gluons for which the scale to has been factored out 
from the components r + , r± (but is still active in the component r_ ) . For the purposes of this work it is enough to 
take for the Lagrangian of the collinear gluons the full QCD Lagrangian and enforce r +1 rj_ <C to when necessary. 



A. Matching QCD to NRQCD-f SCETi 

For definiteness, we shall restrict our analysis to 3 Si states, which decay mainly through two additional gluons. At 
tree level, the electromagnetic current in (JIJ can be matched to the following currents in this EFT pjj 1 

Mx) = e -°"- (rSj^JjJ^Or) + rj 1 *) J(Vso)(x) + r$g>)j$ Pj) (z) + ...) + h.c. (2) 



r (8/5 ) = J%1 S0) (X) = X^BI^X) 

rgif 5 - ^ {ikji* + r&tij ~ ^«') J$Pj)W ^= ~iX iB l VVV(a;) (3) 

where n = (n + , n_, n±) — (1, 0, 0) and = e Q , Atp on p . These effective currents can be identified with the leading order 
in a s of the currents introduced in [JjJ. We use both Latin (1 to 3) and Greek (0 to 3) indices, 73" is a single collinear 



1 One loop matching calculations are already available analytical for the octet currents |2.'I and numerical for the singlet one |24| . 
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gluon field here, and eeQ is the charge of the heavy quark. Note, however, that in order to arrive at one need not 
specify the scaling of collinear fields as M(A 2 , 1, A) but only the cut-offs mentioned above, namely r + , rj_ -C M. Even 
though the P-wave octet piece appears to be 1/m suppressed with respect to the S-wave octet piece, it will eventually 
give rise to contributions of the same order once the bound state effects arc taken into account. This is due to the 
fact that the 3 Si initial state needs a chromomagnetic transition to become an octet So, which is a s suppressed with 
respect to the chromoelectric transition required to become an octet 3 Pj. 
T(z) can then be written as 



where 



and 



T(<A — H^ 1 ' 351 ) rpii'aa'pp' , rr^, 1 So) rpaa' , rr(8, 3 Pj) rpaija'i'j' , , ,\ 

1 \ Z ) - n ii>aa'0P' 1 (l3s 1 ) + V 1 (8, 1 S ) + U aija'i' 3 ' 1 (8.sp J ) "I KV 



rr(l, 3 Si) _ UV r {l, 3 Si) r (l, 3 Sl) 
^ii'aa'pP' ~ '<JL 1 aPin 1 a'P'i'v 



tt^So) _ P^ r (8, 1 S'o) r (8, 1 S'o) 
aa' '/_L L afi L a' v 

tt(S, 3 Pj) ^ ^ r (8, 3 Pj) r (8, 3 P,,) frt 
aijet' i' j' '_L afiij a'ui' j' \ / 
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f d 4 xe 


-iq-x — 2mxo 






T $s )(z) = "* J 


d^xe 


-iq-x — 2mxo 


\ 3 Si\T{J" a i So) 




= -* J 


f d A xe 


-iq-x — 2mxQ 


( K 3 Si\T{J (8 i Pj 





(6) 

In J3J we have not written a crossed term (8, 1 Sq- 3 Pj) since it eventually vanishes at the order we will be calculating. 



B. Matching NRQCD+SCETi to pNRQCD+SCETn 

If we restrict ourselves to z such that M(l — z) < moq, the scale of the binding energy, we can proceed one step 
further in the EFT hierarchy. As discussed in Ref. 21], NRQCD still contains quarks and gluons with energies 
~ ma s , which, in the situation above, can be integrated out. This leads to Potential NRQCD (pNRQCD). On 
the SCET side, the restriction above implies that one may also restrict collinear gluons in the final state to have 
r+,r± <C M VI — z < ma s , as we shall do. The scale My/1 — z, which is still active in SCETi, must then be 
integrated out 0. The integration of this scale produces the dominant contributions from the color singlet currents. 
We have 

( 3 S 1 \T{Jif 3Si) (x) i j{ 1 %(0)}\ 3 S 1 )^ 



2N c S v \ Xl O,x )Sy (0,0,0) (\A.G\Tr{BlBl}{x)Tr{BiB^ }(0)|VAc\ (7) 



The calculation of the vacuum correlator for collinear gluons above has been carried out in [Tlj . and the final result, 
which is obtained by sandwiching between the quarkonium states, reduces to the one put forward in that reference. 
For the color octet currents, the leading contribution arises from a tree level matching of the currents l|2jl. 



,/; 



Or) 



J (8. 3 P 



^/2T^O a P (x,0,x )Bl a (x) 



Bl a (x) 



(8) 



Sy, Oy and O p are the projection of the singlet and octet wave function fields introduced in [2l| to their vector and 
pseudoscalar components, namely S = (S P + S v a l )/V2 and O a = (0% + Opa^ /V2, T F = 1/2 and A^ c = 3 is the 
number of colors. B°^(x) in © are now collinear gluons with r + , rj_ <C M \/l — z < ma s . 



FIG. 1: Color octet contributions. • represents the color octet S-wave current, A represents the color octet P-wave current. 
The notation for the other vertices is that of Ref. El, namely <■:= ; gCF {a ^- a ^ Tr \T b B] and ©: = -Ji — xTr \T b E] . The 

solid line represents the singlet field, the double line represents the octet field and the gluon with a line inside represents a 
collinear gluon. 

C. Calculation in pNRQCD+SCET„ 

We shall then calculate the contributions of the color octet currents in pNRQCD coupled to collinear gluons. They 
are depicted in Fig. ^ For the contribution of the P-wave current, it is enough to have the pNRQCD Lagrangian at 
leading (non-trivial) order in the multipole expansion given in |2l|. For the contribution of the .S-wave current, one 
needs a 1 /m chromomagnetic term given in |25| . 

Let us consider the contribution of the 5*-wave color octet current in some detail. We have from the first diagram 
of Fig. □ 



T$sAz) = -^T(^)^tU^Y aMC f f d 3 x / d 3 x'^„(x')t (x) 



(8,is )W- M '-L y "'Y F \2m) " s ^ u ^ J J UA ™^ A ^ n °W 7(2^)4 k 2 + ie 



d 4 k 



k Q + E n -h + ieJ x , (M (1 - z) - k+)M - + ie \ -k + E, 



1 



(9) 



Ox 



where we have used the Coulomb gauge (both for ultrasoft and collinear gluons). E n < is the binding energy 
(M = 2m + E n ) of the heavy quarkonium, i/>„o( x ) its wave function, and h n the color-octet Hamiltonian at leading 
order, which contains the kinetic term and a repulsive Coulomb potential |21| . cp is the hard matching coefficient 
of the chromomagnetic interaction in NRQCD [2(|, which will eventually be taken to 1. We have also enforced that 
k is ultrasoft by neglecting it in front of M in the collinear gluon propagator. We shall evaluate © in light cone 
coordinates. If we carry out first the integration over fc_, only the pole k- — kj_/A;+ contributes. Then the only 
remaining singularities in the integrand are in the collinear gluon propagator. Hence, the absorptive piece can only 
come from its pole M 2 (l - z) - Mk + = k 2 x . If k + < M(l - z), then k^ ~ M 2 (l - z) which implies fc_ ~ M. This 
contradicts the assumption that k is ultrasoft. Hence, k^ must be expanded in the collinear gluon propagator. We 
then have 



im (i^' So )(*)) = -vT\^)fn ^y aMCf x 



3C 



X J d 3 xj d 3 xV; o ( X ')Vno(x)^: J dk+S (M(l - z) - k+) X 

x r^d^' H vxA - h h f z rt— ) i do) 

Jo VI h - E n + -f- + x I h Q - E n + -f- + x h - E n + -f- + x 



x,x' 



where we have introduced the change of variables |kj_| = ^J2k+x. Restricting ourselves to the ground state (n = 1) 
and using the techniques of reference [2(j we obtain 



Im 



J 
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x / ,l.r(2^ w (0)I s (^±+x)-lU^+x) 



Is(-z- + x) := / dVio(x) 



nit - 



2 ^ V /l °- S l+ 2 +a; /x.O 



7 a s iV c 1/ 2z' /A A 1 - 



7T 2 1-Z' 

where 



1-7— ---i.i--.n-? ( n ) 

1 + z' \ z' z' 1 + z' J J 



^ = ? = A = — ^- (12) 



2 7 to 2 27V C C/ 

(Pi = —m ^ Cf 2^ = — ■*—)• This result can be recast in the factorized form given in 11]. 

( T $s )( z )) = -VT' [ dl + S s (l + )lmJ M (l + - M(l - z)) 



Im 



ImJ„(i+ - M(l - z)) = Tl (Nl - 1) ^S(M(1 -z)-l + ) 

= W ©T * (*M°>'4 + *> - *<T + *>) (13) 

We have thus obtained the S- wave color octet shape function Ss(l+)- Analogously, for the P-wave color octet shape 
functions, we obtain from the second diagram of Fig. ^ 



Im 



(z))^»jr'^y«+fe'5pi(i+) + (nV- ) Sp2{1+)) lmJ M (l + - M(l - z)) 
SpiG+) := J/x ^ w (0)I P ( l -± +x)- + *) 

* p2(;+) := r^jrrbf (^ o(0) " 2iM0)/p( y + + /|>( y + x) i (14) 



where 



fHf+l):= __y^, oW , _^_ t _ v 



2 ' " / x,0 



(2 " x) Aiihr ( 2(1 + z ' )(2 + z ' ] + 3z ' )M +x)+ 2M + A)2+ 

1 f ,, ,w ,9 A(l-2') /A A 1-z'\ 



(15) 



Note that two shape functions are necessary for the P-wave case. The three shape functions above are UV divergent 
and need regularization and renormalization. In order to regulate them at this order it is enough to calculate the 
ultrasoft loop (the integral over k in JJJ) in D-dimensions (leaving the bound state dynamics in 3 space dimensions). 
The UV behavior can be easily obtained by making an expansion of Is and Ip in 1/z' , which is displayed in formulas 
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[3j an d (|A4|) of the Appendix. For the purpose of this section we only need the expansions up to order 1/z' 2 . The 
singular pieces read (D = 4 — 2e) 




^ + )Uo-^^f-+ln^] + ---] (16) 



3tt 2 N c V 

The renormalization is not straightforward. We will assume that suitable operators exists which may absorb the 
1/e poles so that a MS scheme makes sense to define the above expressions, and discuss in the following the origin 
of such operators. In order to understand the scale dependence of (I16|) it is important to notice that it appears 
because the term in the collinear gluon propagator is neglected in (|5J . It should then cancel with an IR divergence 
induced by keeping the term kj_, which implies assuming a size M 2 (l — z) for it, and expanding the ultrasoft scales 
accordingly. We have checked that it does. However, this contribution cannot be computed reliably within pNRQCD 
(neither within NRQCD) because it implies that the k- component of the ultrasoft gluon is of order M, and hence 
it becomes collinear. A reliable calculation involves (at least) two steps within the EFT strategy. The first one is 
the matching calculation of the singlet electromagnetic current at higher orders both in a s and in (k±/M) 2 and 
k + /M . The second is a one loop calculation with collinear gluons involving the higher order singlet currents. Notice, 
before going on, that all divergences (and logarithms) of Ss(l+) and Spi(l+) in (|16fl are sensible to the bound state 
dynamics, whereas those for Sp 2 ('+) are n °t- F° r the former, Fig. |2Jshows the relevant diagrams which contribute 
to the IR behavior we are eventually looking for. We need NNLO in a s , but only LO in the (kj_/M) 2 and k + /M 
expansion. These diagrams are IR finite, but they induce, in the second step, the IR behavior which matches the 
UV of l(TB|) . The second step amounts to integrating out the scale My/1 — z by calculating the loops with collinear 
gluons and expanding smaller scales in the integrand. We have displayed in Fig. 0the two diagrams which provide 
the aforementioned IR divergences. For the latter, the UV behavior of which does not depend on the bound state 
dynamics, we need the matching at LO in a s (last diagram in Fig. [21 but NLO in k + /M and (kj_/M) 2 . We have 
checked that the coefficient of the logarithm coincides with that of the (1 — z) log(l — z) term in the QCD calculation 
as it should. 

The above means that the scale dependence of the leading order contributions of the color-octet currents is of the 
same order as the NNLO contributions in a s of the color-singlet current, a calculation which is not available. One 
might, alternatively, attempt to resum logs and use the NLO calculation [24| as the boundary condition. This log 
resummation is non-trivial. One must take into account the correlation of scales inherent to the non-relativistic system 
|27|. which in the framework of pNRQCD has been implemented in [2sH29j ]. and combine it with the resummation of 
Sudakov logs in the framework of SCET (see also H^)- Correlations within the various scales of SCET 

may start playing a role here as well [3lJ . In any case, it should be clear that by only resumming Sudakov logs, as it 
has been done so far Q , one does not resum all the logs arising in the color octet contributions of heavy quarkonium, 
at least in the weak coupling regime. 



III. APPLICATION TO THE T(1S) 

We apply here the results of section[n]to the Y(1S). There is good evidence that the T(1S) state can be understood 
as a weak coupling (positronium like) bound state p^L flfi! UtIIlrI ] . Hence, ignoring O ( Apnn) in the shape functions, 
as we did in section [HI should be a reasonable approximation. In the analysis of the effects of the octet shape 
functions were set to zero, so we expect to improve on their results. We plot in Fig. 0]the CLEO data in the end 
point region |l2j |. the curve obtained in hjl (dashed line) and our curves (solid and dot-dashed lines). Our curves 
are obtained by adding to the results of 11], which consist of the LL resummation of the singlet contributions only, 
our (MS) results for the color octet contributions (without LL resummation 0) and setting the scale dependence to 
fi = My/l — z (solid line) and to \x — 2 ±1 MVl — z (dot-dashed lines). The first choice is the most reasonable one 
according to the discussion in the previous section and the last ones are displayed in order to get the flavor of the 
systematic errors. We have used the following values for the masses and a s in our plots: mi, = 4.8 GeV, My = 9.46 
GeV, a s (fih) = 0.216, a s (/i s ) = 0.32 and a s (fi u ) = 0.65. [ih ~ m stands for the hard scale and is to be used for the 
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FIG. 2: Relevant diagrams in the matching calculation QCD -> pNRQCD+SCET. 





FIG. 3: Diagrams which induce an IR scale dependence which cancels against the UV one of the octet shape functions. 



a s arising from (3). /i s ~ ma s stands for the soft scale and is to be used for the a s participating in the bound state 
dynamics (in E±, Is, Ip and ?/>iq(0)). /j, u ~ ma 2 stands for the ultrasoft scale and is to be used for the a s arising 
from the coupling of ultrasoft gluons. It turns out that the color octet contribution is numerically enhanced and 
dominates over the color singlet one in the whole end-point region. In order to compare with the experimental curve, 
the theoretical result must be convoluted with the experimental efficiency |T^ and the overall normalization must be 
taken as a free parameter. By adjusting our curves to data around z ~ 0.7, we obtain an almost perfect agreement 
in the whole end-point region (z € [0.7, 1]) for /i = M\J\ — z (solid line). A complete analysis, including systematic 
errors, is beyond the scope of this paper. It would require either a NNLO matching or a NLO one |24| with NLL 
resummation of the singlet current. In addition one should estimate what the leading non-perturbative effects are. 
In any case, it should be clear from our results that the introduction of a gluon mass [32l | is not necessary for the 
description of the experimental data on the photon spectrum in the end-point region of T(IS') radiative decays. 



IV. DISCUSSION 



We would like to make a few remarks which stem from the details of our calculation. In the existing formulations of 
SCET, suitable scaling properties are assigned to the various modes. However the standard assignments are violated 
in our case. The collinear gluons in Fig. scale as m(A 2 ,l,A 2 ) rather than m(A 2 ,l,A) (A = \f\ — z), which is the 
assigned scaling in 0, 0, ■ This indicates that SCET should better be discussed in terms of UV (and IR) cut-offs for 
the relevant modes, rather than scaling properties. This becomes particularly clear when we analyze the scaling of the 
ultrasoft gluon in the same diagrams. If ma 2 ~ M(l — z), it scales like m(A 2 , A 2 , A 2 ) (A = a s ), which coincides with 
the standard scaling rules. If, however, ma\ ~ M(l — z), it scales like m(A 4 , A 2 , A 3 ), the typical scaling of the recently 
discovered soft-collinear modes Either scaling is properly described by the ultrasoft gluons of pNRQCD, since 
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fi = My/ (1 — z) (natural choice) and the dot-dashed lines are obtained by setting fi — 2My/l — z and /j, — 2 Ms/l — z. 



they are defined as the ones having all four momenta much smaller than the soft scale ~ ma s (UV cut-off), which 
is fulfilled in both cases. However, if one insisted to assign to the ultrasoft gluons a momentum scaling m(A 2 , A 2 , A 2 ) 
and not smaller (i.e. one is introducing an IR cut-off for them), then the situation ma* ~ M (1 — z) would require 
the introduction of new (ultra)soft-collinear modes scaling like m(A 4 , A 2 , A 3 ) with an UV cut-off ~ mA 2 for the + and 
_L components. Whether it is convenient or not to make such splitting is a matter of debate 0, Q • 

The case ma s ~ M(l — z) has not been discussed. For the color-octet contributions, it requires a calculation in 
NRQCD, since if one attempts at doing it from Fig. ^ one immediately realizes that the four momentum of the 
ultrasoft gluon is ~ ma s , a region where pNRQCD is not applicable. Hence at NRQCD+SCETi level (Section [OJ 
one should calculate a set of diagrams involving a collinear and a soft gluon. The leading order contribution comes 
from the 5-wave current only and we have seen it to vanish. 

We have refrained ourselves to put forward a Lagrangian for SCET which also holds for heavy quarkonium systems 
because several issues, like the remarks made above, should be better understood. Clearly, as we have shown in 
this paper, one cannot simply take over the SCET for heavy-light systems and apply it to heavy quarkonium (for 
instance, one misses logs which depend on the binding effects). Our analysis also indicates that it may be convenient 
to rephrase SCET in terms of cut-offs rather than in terms of scaling properties of the various modes as it has been 
done so far. Then, one would have to account properly for both the suitable cut-offs of pNRQCD and those of SCET. 

V. CONCLUSIONS 

Apart from making a few remarks, which we hope will be useful for an eventual construction of a SCET Lagrangian 
adapted to heavy quarkonium systems, we have calculated the S- and P-wave octet shape functions in the weak 
coupling re gime . We have also discussed their scale dependence. The addition of these contributions to the ones 
obtained in makes the agreement with data for the end point photon spectrum of inclusive T(1S) decays almost 
perfect. As a by-product the NRQCD matrix elements <T(lS , )|0 8 ( 1 S , o)|T(lS')) and (T(1S)\0 8 ( 3 Pj)\T(1S)) have 
also been calculated, in the weak coupling regime, (see Appendix^) and their scale dependence discussed. 



FIG. 5: Diagrams which require a Vi( 3 Si) operator for renormalization. The solid circle stands for either the Os( 1 So) or 
Og( 3 Pj) operator, the crossed box for either the chromomagnetic ) or chromoelectric (O ) interaction in Fig. 1, the empty 
box for the octet Coulomb potential, and the thin solid lines for free QQ propagators. 





FIG. 6: Diagrams which require a Oi( 3 5i) operator for renormalization. Symbols are as in Fig. 
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APPENDIX A: CALCULATION OF T(1S) NRQCD COLOR OCTET MATRIX ELEMENTS 

The calculation in Section Til CI can be easily taken over to provide a calculation of (T(1S')| Og( 1 'S'o)|T(l»S')) an d 
(T (1S)\Os( 3 Pj)\T (IS)) assuming that ma 2 ^> Aqcd is a reasonable approximation for this system. Indeed, we only 
have to drop the delta function (which requires a further integration over k+) and arrange for the suitable factors in 
(jT3|l and l|Tt|). We obtain 

<T(15)|O 8 ( 1 5 )|T(15)) = -2T 2 F {N 2 C - 1) / dk + S s (k+) 

Jo 

(T(lS)\O s ( 3 Pj)\T(lS)) = _ 4(2J+ 1)^(^-1) J™ dk+Spi[k+) (M) 

where we have used 

/ dk+S P2 (k + ) = - / dk + S P1 (k+) (A2) 
Jo J Jo 
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The expressions above contain UV divergences which may be regulated in the same way as in Section III CI namely 
by calculating the ultrasoft loop in D dimensions. These divergences can be traced back to the diagrams in Fig. [S] 
and Fig. El Indeed, if we expand Is and Ip for z' large, we obtain 



7 a sN c 

7T 2 



I + i f (-l + 2Aln2) 



1 



1-2A 



A 2 7T 2 



+ -j -1 + A(21n2 + 1) + A 2 (-41n2) + -C(3)A 3 + 0(-g) 



(A3) 



2z' V 4 
1 1 



1 - A + —(-6 + tt z )\ z ] — + - (-5 + A + A 2 (2 - 8 In. 2) + 8Aln2+ 



+A 3 (■ 



-41n2 + 3C(3))) -L + 0(-L)} 



(A4) 



It is easy to see that only powers of 1/z' up to order four may give rise to divergences. Moreover, each power of 1/z' 
corresponds to one Coulomb exchange. Taking into account the result of the following integral, 



/ dk+ / dx{2k + x)- £ — = X- 
Jo Jo z \m 



2-2e 



r 2 (i-e) 
r(f) 



/a 
V2 



2e- 2 



(A5) 



we see that only the 1/z' 2 and 1/z' 4 terms produce divergences. The former correspond to diagrams in Fig. 0and the 
latter to Fig. H3 which can be renormalized by the operators Vi( 3 Si) and Oi( 3 Si) respectively. It is again important to 
notice that these divergences are a combined effect of the ultrasoft loop and quantum mechanics perturbation theory 
(potential loops ) and hence it may not be clear at first sight if they must be understood as ultrasoft (producing 
log/z u in the notation of Refs. [28l l^ - ) or potential (producing log/i p in the notation of Ref. |H,|2^|). In any case, 
the logarithms they produce depend on the regularization and renormalization scheme used for both ultrasoft and 
potential loops. Notice that the scheme we use in this work is not the standard one in pNRQCD 17, 28, 3j| . In the 
standard scheme the ultrasoft divergences (anomalous dimensions) are identified by dimensionally regulating both 
ultrasoft and potential loops and subsequently taking D — > 4 in the ultrasoft loop divergences only. If we did this 
in the present calculation we would obtain no ultrasoft divergence. Hence, in the standard scheme there would be 
contributions to the potential anomalous dimensions only. The singular pieces in our scheme are displayed below 



(T(15)|O 8 ( 1 5 )|T(l^))|^ 



e \ [ira 



-2s 



C F N c a s (ix v 



(C/a s (^ s )) 4 \ 

7T 



2+ 



-A [-7 -4 log 2] + A 2 



4 + 81og2 + 41og 2 2+ — 
3 



+ A 3 



-41og 2 2-^-^C(3) 



(T(1S)\Os( 3 Pj)\T(1S))\^ o,-(2J+1)- S ± Cf aM(C f aMfp, 



x (2 — A) —4 + +A 



47 

^ + 51og2 



+ A 2 



5 2tt 2 8, „ 8, n 
log 2 log 2 2 

6 9 3 8 3 s 



+A 3 



(A6) 
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